Two-level atom at finite temperature 
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Properties of a two-level atom coupled to the quantized electromagnetic field at finite temperature 
are determined. The analysis is based on a new method (inspired by QED) of describing qubits, 
developed previously at zero temperature In this paper, we make a generalization to finite 
temperature by introducing the Matsubara formalism and the temperature propagators. We analyze 
the spectral properties of different types of propagators and we derive a direct connection between 
the temperature propagators and the real time propagators. To show the effectiveness of this 
method, we calculate the temperature dependence of the polarizability of a two-level atom in the 
lowest order of perturbation theory and we predict an unexpected sharpening of the resonance. The 
whole discussion is carried out without the rotating wave approximation. 
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I. INTRODUCTION 

The main purpose of this paper is to show that a 
new method of describing qubits within the quantum 
field theory formalism developed by us in [l| can be also 
used to describe qubits at finite temperature. The main 
idea comes from the observation made long time ago 
by Matsubara that there exists a close analogy be- 
tween Feynman propagators and the so called tempera- 
ture propagators. In his paper, Matsubara described gen- 
eral rules of perturbation techniques at non-zero temper- 
ature based on Feynman diagrams and used them with 
great success to describe the system made of electrons 
and positrons. The Matsubara formalism was further 
developed by many authors [1, 0, d, Q and it has quickly 
become one of the main tools of modern statistical me- 
chanics 0, H, • In the usual treatment the methods 
of quantum field theory are exploited to describe systems 
in which the number of particles is not conserved, i.e. in 
the grand canonical assemble. In this paper we will adopt 
this formalism to analyze the temperature properties of 
one qubit realized by precisely one electron in atomic 
states. It is not obvious that it is indeed possible. There- 
fore, at the beginning we present a short introduction to 
thermal field theory and we give a clear prescription how 
to proceed. 

The influence of thermal fluctuations on the properties 
of qubits was partially discussed before. First, a very 
general treatment was based on dissipation phenomena 
in the damped cavity model [ll|, [l^, Il^l ■ The problem 
of the influence of the temperature on the induced in- 
teractions between qubits via the electromagnetic field 
and on their entanglement is recently considered in the 
context of the so called sudden death of entanglement" 
phenomenon [ij, [l^ [l^ . Very often this analysis is car- 
ried out in the Jaynes-Cummings model in the rotating 



wave approximation (RWA) [T3, [H, [H, US] • There exists 
also some numerical predictions on the behavior of qubits 
at finite temperatures realized in quantum dots [22]. So 
far a general problem of the dynamical properties of a 
single qubit interacting with the quantized electromag- 
netic field (without RWA) at finite temperature has not 
been discussed since all current techniques did not give 
such possibilities. In this paper we want to show that 
methods ingrained in quantum field theory could be the 
first step to solve this problem. 

To make all our argumentation clear, we restrict our- 
selves to the problem of a two-level atom and its polar- 
izability (response tensor), since a two- level atom is the 
simplest model of the interaction between light and mat- 
ter where purely mathematical problems are as small as 
possible Moreover, the properties of the linear po- 
larizability of a two-level atom were discussed in detail 
at zero-temperature regime by many authors and many 
discussions on it took place [ll, H [21, [H, [13, IH, SI, ■ 
This discussions have never went beyond zero tempera- 
ture case, since it seemed to be a very hard task. In this 
paper we want to show that our new method of describ- 
ing qubits in terms of field operators is a key to manage 
this problem. To do it correctly we will not make RWA 
at any point of the analysis j28j . 



II. PHYSICAL SITUATION 



A. The Hamiltonian of the system 



The starting point of our analysis is the well known 
Hamiltonian of the two-level atom interacting with the 
quantized electromagnetic field 
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"H = Hq + Hi 



(1) 
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where 



C. Statistical operator 



Ho — mcTz + / dkka\k)a{k), 
Jq 

/>oo 

Hi ^(7^ / dfc.g(A:)$(fc). 



(2a) 
(2b) 



The Hamiltonian Tip describes the properties of noninter- 
acting subsystems - two-level atom with the energy gap 
equal 2m and the free quantized electromagnetic field. 
The operators a^{k) and a(k) create and annihilate pho- 
tons in the appropriate modes of the electromagnetic field 
and together they form a quantum scalar field 



$(fc) 



a(fc) -f at(fc) 



2k 



(3) 



The Hamiltonian Tii is responsible for the interaction 
between the photons and the two-level atom. The form- 
factor g{k) is simply a coupling coefhcient measuring the 
strength of a dipole transition in the atom. Note that 
we do not make a rotating wave approximation in the 
coupling Hamiltonian (j2bp . As was shown previously in 
(28t the complete interaction term is necessary to obtain 
a correct expression for the polarizability of a two-level 
atom. 



B. Two- level atom in thermal equilibrium 

In our treatment of a qubit at non-zero temperature 
we will make two physical assumptions. First, that our 
system described by the Hamiltonian ([T]) is in thermal 
equilibrium with the heat reservoir (thermostat) at tem- 
perature T. Second, that the temperature is not so high 
as to induce the creation of real electron-positron pairs. 
It means that the number of fermions in our system is 
conserved. Notice however, that we do not assume that 
the number of photons in the system is conserved. 

Description of the statistical properties of a quantum 
system in thermodynamic equilibrium is realized in terms 
of the statistical assembles. In our model an appropriate 
assemble is defined by the following control parameters: 

• the number of fermions (electrons) equals 1, 

• the temperature T of the thermostat, 

• the chemical potential fi of the photon reservoir. 

It is worth to notice that because of the properties of the 
electromagnetic field, the chemical potential fi is always 
equal to 0. It means that all statistical properties of the 
photons can be expressed as functions of the temperature 
T and one extensive parameter (for example the volume 
of the cavity). Nevertheless, to emphasize the fact that /i 
is one of our control! parameter we will write it explicitly 
to the end of this section. 



From the rules of quantum statistical physics (cf. , for 
example, ^) we know that any system in equilibrium 
with the thermostat and the reservoir of particles is in 
the quantum state represented by the following density 
matrix (in this context also called statistical operator) 



1 



-I3K 



Z = tr 



(4a) 



As always = l/kT and Z is partition function. In our 
case the statistical Hamiltonian K, is given by 



(4b) 



where the operator A/" represents the number of photons 
in the system. Expectation value of any operator O in 
such a quantum state of the system is given by 



((O)) =tr[pO] 



tr e" 



o] 



tr[e 



-/3/cl 



(5) 



For future needs let us now introduce an additional 
system composed of a noninteracting two-level atom and 
the electromagnetic field. We define the free statistical 
Hamiltonian /Co and the free statistical operator po of 
such a system as follows 



1 

Po = — e 
2o 



Zo = tr 



(6a) 
(6b) 



The statistical operator po describes a quantum state of 
a system composed of a free two-level atom and the free 
electromagnetic field which together are in equilibrium 
with the thermostat and the reservoir of photons. The 
connection between statistical Hamiltonians in the inter- 
acting and the noninteracting case is following 



/C = /Co + Til 



(7) 



and the expectation value of the operator O in the state 
Po is given by 



((0»o-tr[poO] 



tr[e 



-13K.0 



(8) 



D. Linear polarizability of the atom 

The linear polarizability of a two-level atom is a phys- 
ical measure of the reaction of the system to small exter- 
nal electromagnetic perturbations. When the state of the 
system is described by the statistical operator (j4a| , then 
the polarizability is defined as the following two-point 
retarded correlation function [33, ^Sl 



a{t,t') = -iA6{t-t') 



e^«V,e-'^*,e^«*a,e- 



mt' 



(9) 
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For a single two-level atom, the constant A is equal to 
(P /2h where d is a dipole-moment of the atomic transi- 
tion. It is easy to find that a{t, t') depends only on the 
difference of its arguments and therefore it has a very 
simple form in the frequency domain 

/oo 
dte*"(*-*')a(i,i')- (10) 

In this paper we give a precise and clear prescription 
how to find this quantity as a perturbation series using 
the well known methods of thermal quantum field theory 

mMMM- 



III. SECOND QUANTIZATION 

The basic ingredients of quantum field theory are ob- 
viously the field operators. Therefore, before we start to 
analyze our system in this language, we have to intro- 
duce, besides the electromagnetic field operator ([3]), the 
fermionic field operator representing the states of an elec- 
tron in the atom. As usual, field operators can be intro- 
duced after performing second-quantization of electronic 
states. In this model this procedure is very easy. We just 
enlarge the two dimensional Hilbert space spanned by 
vectors |g) and |e) representing ground state and exited 
state of the electron in the atom to the four dimensional 
one spanned by four state vectors: two previously defined 
|g) and |e), the zero-electron state |N), and the antisym- 
metric two-electron state |B) representing two electrons 
- one in each state. 

Such a procedure enables one to introduce creation and 
annihilation operators of the electron. The operators ipi 
and f/'J annihilate and create electron in the ground state 

of the two-level atom, and ipi and ip] do so in the exited 
state. These operators form together the fermion field 
operator and its hermitian conjugate [l| 

(11) 

Then one can easily find the second-quantized form of the 
new Hamiltonian acting in the extended Hilbert space: 

H = Hq+Hi (12a) 

where 

/>oo 

n^^m^^a^^+l dkka){k)a{k), (12b) 

^0 

/•oo 

= ^-to-^r * / dfc.g(fc)$(fc). (12c) 
Jo 

The origin, and the properties of the Hamiltonian ([72]) 
and of the field operators ^ and are discussed in de- 
tail in our previous paper Notice that the system 
described by the Hamiltonian TL may have states with 



different number of fermions. However, the number of 
fermions operator commutes with 7i. 

In the Heisenberg picture, the time evolution of the 
system is contained in the dynamics of the field operators. 
They evolve according to the following equations [l|: 

/>CXD 

(ift - ma^) *(t) = / Akg{k)^{k,t)<j,-^[t), (13a) 
Jo 

(9,2 + e) $(fc, t) = -g{k) ^>\t)cj,^{t). (13b) 

Now we introduce a new statistical operator p repre- 
senting the thermal equilibrium of the new quantum me- 
chanical system 

p = ^c-'^'^, Z = tr [e"'^'^] (14a) 

where the statistical Hamiltonian /C is given by 

K, = n- pM. (14b) 

The expectation value of any operator O in the quantum 
state p is given by 

Notice, that the expectation values ^ and ([T5|) . even 
for the same physical observable O, are different since 
the statistical Hamiltonians K, and /C act in different 
Hilbert spaces. Nevertheless, for those qubit operators 
d , . . . , On that are represented by traceless 2x2 ma- 
trices (for example the free qubit Hamiltonian moCTz, the 
interaction Hamiltonian (j2b[) . or any other linear combi- 
nation of the Pauli matrices) one can verify that there 
exists the following direct connection between the ex- 
pectation value of their product ((Oi • • • On)) and ex- 
pectation value of their second-quantized counterparts 

{{Oi---On)) = |(*tOi*---*ta„^'). (16) 

It should be emphasized that this relation is a central the- 
orem of this paper. It expresses nontrivial but a very sim- 
ple connection between the expectation values of physical 
quantities of the system containing exactly one fermion 
and the expectation values of their non-physical second- 
quantized counterpart. Therefore, it is a peculiar link 
between the interesting physical system described by the 
Hamiltonian ([1]) and a convenient mathematical model 
described by the Hamiltonian (H^). Without this con- 
nection further analysis would not be possible. 

Relation follows from the observation that for 

traceless matrices the bilinear combination of field oper- 
ators ^''l'Oj\E' gives zero when acting on the two-electron 
or the zero-electron states. Moreover, in the qubit (one- 
electron) subspace it acts as the qubit operator Oi. Ob- 
viously, this argument is also valid for evolving qubit 



4 



-mt 



and their counterparts 



operators 0^{t) = e^'^^de 

This all means that a linear polarizability of the two- 
level atom defined by traceless matrices in © can be 
reproduced from its second-quantized counterpart 

(17) 

by using the property p6)) . 

Note that function a{t, t') has no direct physical sense, 
since the statistical operator p involves also non physical 
states. Therefore we stress that we treat a{t,t') only as 
a very convenient tool to find a physical quantity. 



IV. MATSUBARA FORMALISM 

The starting point of the formulation of quantum field 
theory at finite temperatures is to introduce the temper- 
ature Green functions (or temperature correlation func- 
tions). They are defined in an analogous way as the Feyn- 
man Green functions at zero temperature. The difference 
lies in the evolution parameter - instead of the ordinary 
time variable t we use the imaginary time r — it. As was 
first noticed by Matsubara in 1955 after this simple 
replacing one can use almost all known machinery of the 
quantum field theory. In analogous to the ordinary time 
case all analysis is carried out in two quantum mechanical 
pictures - Heisenberg-like picture and Dirac-like picture. 

In the Matsubara-Heisenberg picture the whole evolu- 
tion of the system is contained in the evolution of the 
operators. Any operator O evolves according to the fol- 
lowing rule 



0{t) = c^''Oe-'^\ 



(18) 



In particular, the field operators Q and (jlip evolve in 
the following way 



$(fc,r) = e'^"'$(fc) 



*(t) = c'^^^'c 
*t(r) = c'^^'^'te 



-Kt 



(19a) 
(19b) 
(19c) 



In contrast, in the Matsubara-Dirac picture the oper- 
ators evolve according to the free statistical Hamiltonian 
/Co 

0{t) = e'^"'' O 6"^°^ (20) 

One can easily check that the field operators ^ and ITTI) 
evolve in the following way (from now we will put /i = 
everywhere) 



V.t(r) 

0(fc,T) 



mr 
7nr 



( 



e™^V^e-"^), 
a(fc)e-'=^ + at(fc)e'=^ 



2k 



(21a) 
(21b) 
(21c) 



The interaction Hamiltonian (|12cp in this picture has the 
form 



-KaT 



— w.t 



i;^{T)a^,p{T) I dk g{k)(t>{k,T). (22) 





A. Temperature correlation functions 

Temperature correlation functions are the central ob- 
jects of the quantum field theory at finite temperatures. 
In analogy to the Feynman ones, they do not have any di- 
rect physical interpretation, but they are the basic ingre- 
dients of the perturbation computations based on Feyn- 
man diagrams. 

We call the n-point temperature correlation function a 
complex function of n parameters defined as follows (To 
underline the difference between imaginary and ordinary 
time functions we use a ,, tilde"): 



g(ri,...,r„)--(T,0(ri)...0(T„)). 



(23) 



where the operator T,- is the time-ordering operator. 
This operator respects the fermionic nature of operators, 
and therefore for any odd permutation of electron opera- 
tors one should remember about changing the sign of the 
expression. It will be clarified later that in practice it is 
sufficient to restrict each imaginary-time variable only 
to the range 



(24) 



0<n<(3. 



Special case of temperature correlation functions 
are temperature propagators of electromagnetic and 
fermion fields. They are defined as follows 

5(1/3 (ti - T2) = Saf}{Ti,T2) 

= -(T,*„(ri)*]j(T2)), (25a) 

= -(T,$(fc,Ti)$(A:',T2)). (25b) 

The temperature propagators are functions of the dif- 
ference of their arguments only. Therefore, according to 
the restriction (IM)) they are defined on the finite range of 
length 2/3. Moreover, they fulfill the symmetric boundary 
condition Q{—0) = GiP)- Therefore one can represent 
them as the following Fourier series 0, H, [1, '■ 



1 °° 



' S(w„), 



(26a) 



where 



2^ bosons. 



2(n+l)7r 

P 



fermions. 



(26b) 
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The Fourier coefficients Q (u;„) can be easily reconstructed 
from the inverse Fourier transform 



gK)= / dre^""^e(T). 
Jo 



(27) 



Since all properties of a two-level atom will be deter- 
mined via the free temperature propagators we will now 
calculate them in detail. Similarly as it was done at zero 
temperature [l|, free propagators are defined in terms of 
the field operators in the Matsubara-Dirac picture: 

So,f3(Tl,T2) = -(T,Va(ri)^^(T2))o, (28a) 

D{k,k',Ti,T2) = -{Tr<t>{k,Ti)ct>{k',T2))Q. (28b) 

To find the free temperature electron propagator one 
has to use the averages of the bilinear combinations of 
the fermionic field operators over the statistical assem- 
ble. Then by simple calculations one finds the following 
representation of the electron propagator: 

g, X _ _ ^(^1 - ^2) _ 0{ri-T2) \ 

(29) 

Consequently the Fourier coefficient of this propagator 
reads 

1 , (2n-|-l)7r 



where 



(30) 

Similarly, by using appropriate averages of bosonic 
field operators, one can get the free propagator of the 
electromagnetic field. It is not hard task to find the 
free temperature photon propagator and to show that 
its Fourier component reads 



D[k, k ,LUn) ^ ^ 



where 



2n7r 
(31) 



B. Perturbation series and Feynman rules 

The methods of quantum field theory at finite temper- 
ature are based on theorems which are formally very sim- 
ilar to the standard theorems used at zero temperature. 
This observation gives us a possibility to adapt automat- 
ically all well known machinery to this more complicated 
physical situation. 

The main relation which gives us a chance to start 
perturbation calculations is a generalized Gell-Mann and 
Low theorem. It gives us explicitly a connection between 
the temperature correlation functions in the Matsubara- 
Heisenberg picture and the corresponding correlation 
functions in the Matsubara-Dirac picture. This funda- 
mental relation is [7, 9] 

(T,Oi(Ti)...0„(r„)) 

_ (T,Oi(ri)...0„(T„)e-/o' 



(T,e-/o'^<i-«'M)o 



(32) 



One can easily check that the formal structure of this re- 
lation corresponds to the structure of the corresponding 
equation at zero temperature [H . Since the integration on 
the right hand side is done only over the range [0, /3], now 
it is clear why we restricted in (j24[) the imaginary-time 
parameters to this range. Formal expansion of the right 
hand side in the powers of interaction Hamiltonian leads 
us directly to the concept of Feynman diagrams which are 
just simple graphical representations of the appropriate 
elements of the perturbation series. As usual any temper- 
ature correlation function in the Matsubara-Heisenberg 
picture can be expressed as a infinite sum of connected 
Feynman diagrams. 

The finite-temperature Feynman rules for studied 
qubit system are very simple. At first one should draw 
Feynman diagrams according to the same rules as were 
used at zero temperature [1]. Then one assigns them cor- 
responding mathematical expressions. It can be done in 
the frequency domain by applying following rules: 

• Each photon line between any two vertices repre- 
sents the free temperature photon propagator 

I.' ^ 5{k-k') 

where ujn — 

• Each fermion line between any two vertices repre- 
sents the free temperature fermion propagator 



-S{uJn) 



where w„ = -^^ — ■p-'— . 



-1 



• Each vertex represents the expression 

~V{k) - -g{k)a,. 

• In each vertex the frequency conservation law 
holds. It is a consequence of the imaginary time 
translation invariance of the system. This rule is a 
counterpart of the energy conservation law at zero 
temperature. 

• The fermion propagators are 2x2 matrices and 
therefore they should be multiplied in the order 
given by the directions of the lines. 

• Each closed fermion loop leads to the trace oper- 
ation of the matrices and gives a factor —1. This 
rule is a simple consequence of the fermionic nature 
of the electron filed operators. 

In the end one should integrate such an expression over 
all internal momenta k and take a sum over all internal 
frequencies uj„. Moreover each sum must be divided by 
/3. 
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V. TEMPERATURE PHOTON PROPAGATOR 

Now we will find the temperature photon propagator in 
the lowest order of perturbation theory. Since there exists 
a formal connection between all Feynman rules at finite 
and zero temperatures, we can automatically write down 
the fundamental relation between photon propagator in 
the Matsubara-Heisenberg and in the Matsubara-Dirac 
pictures. This relation reads [l| 



dki / dfc2 D{k, k , u;„)n(fci, fc2, w„)I?(fc2, k , w„). 
Jo 

(33a) 



We introduced here the so called temperature self- 
energy function 11 which is built from all irreducible con- 
nected Feynman diagrams with two external photon lines 



Il{kuk2,U;r^)=0 = Q) + (h) 



(33b) 

One can transform the relation (j33p to a more convenient 
form 



V[k,k',uj^) 



5{k-k') 



fc2 



gik') 



where we introduced a new function 



h{uJn) = - 



dfc- 



fc2' 



fc'2' 

(34a) 



(34b) 



The summation oven n' can be easily done and the details 
of these calculations are presented in Appendix B. We 
rewrite a result (IB5I) derived there as follows 



4m 



■ tanh 



/9m 



(37) 



4m2 + cj2 \^ 2 

From the cq. (I35p one can easily find the temperature 
transition matrix 



4m 



(4m2 + ujI) coth j - Amh{uJr^) 



(38) 



It is worth to notice that taking the zero-temperature 
limit (/3 — > oo) at this point of our analysis has no phys- 
ical sense. This fact will be explained later. 



VI. REAL TIME PROPAGATORS 



Let us now return to the problem of the linear po- 
larizability of the atom and also introduce the so called 
retarded photon propagator. For the operators in the 
Heisenberg picture this propagator is defined as follows: 



V^{k,k',t,t')^-ie{t-t')mt),<i>{t')]) (39a) 



We also define the temperature transition matrix T 
1 1 



T(tJ„) 



(35) 

Now we are ready to find the lowest order correction 
to the temperature Feynman photon propagator. The 
temperature self-energy function of the photon in this 
order is represented by the following Feynman diagram 



.p(2)(. 




^Tr{ crabs' (a;„' + UJn)cr^S{Ldn')] ■ 



(36) 



and for the free photon field (in the interaction picture) 
it reads 



D^{k, k', t, t') = -ie{t - t'){[cbit), 0(i')])o. (39b) 



There exists a general connection between the retarded 
photon propagator defined above and the polarizability of 
the atom given by ([T7]). This connection follows from the 
equation of motion (jl3bp , a definition of the polarizabil- 
ity ^T7\ . and the properties of the free retarded photon 
propagator (|39b[) . The leading argument is identical as 
at zero temperature and it was given by us in [J] . There- 
fore we will not repeat it here. This connection has a 
form 
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1 r°° r 

2?R(fc, k', ko) = Dj,{k, k\ ko)-^J dfci g{ki) J 



dk2 g{k2) Dfi{k, ki,k(i)a[ko)D-s.{k2,k' , k^) 



(40) 



where 'Dj^{k , k' , ko) and D^{k,k\kQ) are the Fourier 
transforms of the retarded photon propagators defined 
by eq. ^ 

/OO 
die*'=°*I?R(fc,A:',t,0), (41a) 
-OO 

/OO 
dte'''°^ D^{k,k',t,Q). (41b) 
-OO 

From the eq. (PD|) it clearly follows that we need an 
appropriate retarded photon propagator to find the po- 
larizability of a two-level atom (or more generally - the 
linear response of the system). 



A. Spectral representation of propagators 

Now we will demonstrate that there exists an exact 
mathematical connection between the real time and the 
temperature propagators. For future benefits we restrict 
ourself just to the retarded propagator of the electro- 
magnetic field. Other propagators, like the advanced 
propagator or the chronologicaly ordered propagator, can 
be also reconstructed from the temperature propagators 
0, H, but in the case of the linear polarizability of 

the atom they are not interesting. 

From the equation (|27p one can find Fourier transform 
of the temperature photon propagator defined by (j25bp 



J 



V{k,k',LUn)^ dre*'^"^X>(/fc,fc',r,0) 



f0 1 

/ dre*'^"^- V(r^|e-'^'^e'^^$(fc)e-'^^|m)(TO|4'(/^')|f^) 
Jo ^ 



= J2 dre('""+^"-^'")^- 



(n|$(fc)|r7i)(m|$(fc')|n) 



(42) 



r 



In this sequence of equalities we have used the set {\n)} 
of the eigenvectors of the statistical Hamiltonian K. with 
the eigenvalues Kn- After performing an integration over 
T, one finds 



V{k,k',UJn) 



°°dM»M), (43) 



where the temperature spectral matrix of the propagator 
is given by 



m{M, k, k') 



-I3K„ 



X (n|<E>(yfc)|m)(m|$(fc')|n). (44) 

In an analogous way one can show that retarded prop- 
agator (|41ap has the following spectral representation 



Vj,{k,k\t,0) = -i0(t)[($(fc,t)$(fc')> - ($(fc')$(fc,i)) 



2n 



M - 



(45) 



It means that the Fourier transform of the retarded pho- 
ton propagator reads 



V^{k,k\ko)= f dM^ 



miM, k, k') 



M- 



(46) 



B. The analytic continuation 



The spectral representation of the propagators is very 
useful to find the connection between different types of 
propagators in the frequency domain. First notice that 
the retarded propagator treated as a function of the 
continuous complex variable ko is analytic in the upper 
half of the complex plane. It means that if one defines 
the function F{z) of the complex variable z, 



F{z) 



AM 



m{M, k,k') 
z-M 



(47) 



then this function in the upper halfplane would be ex- 
actly equal to the retarded photon propagator ([46]) . Ad- 
ditionally, at points z„ = iuin this function is equal to the 
Fourier components of the temperature photon propaga- 
tor (|43p . Therefore the temperature photon propagator 
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for positive aj„ can be easily reconstructed directly from 
the retarded propagator in the following way 

■D{k,k',ujn)^Vn{k,k',iiUn). (48) 

The fourier components for the negative frequencies can 
be reconstructed as well because the following relation 
holds 

V{k,k',-LOn)=V*{k,k',LUn)- (49) 

The inverse problem, how to construct the retarded 
propagator from the temperature propagator, is more 
complicated. The reason is that the Fourier components 
of the temperature propagator are defined only on a dis- 
crete set of frequencies iVn- Therefore, there exist many 
analytical continuations of the temperature propagator 
to the whole complex plane and one needs to add other 
conditions to choose an appropriate continuation. This 
general problem was discussed in the past 0,B[1,[33| and 
now it is known how to manage it. In our case the re- 
tarded propagator can be constructed with the following 



rule 

Vnik,k',ko)^ F{~iko + e), for ko > 0. (50) 

To find the retarded propagator for the negative fcg one 
has to use the crossing relation given in our previous 
paper 



VII. POLARIZABILITY OF THE TWO-LEVEL 
ATOM 

In the previous section we have found a connection 
between the temperature photon propagator and the re- 
tarded photon propagator. It was very important step 
because it gives us a possibility to find the temperature 
properties of the linear polarizability of the atom. To 
make it obvious let us recall once more a fundamental 
formulas for different types of the photon propagators in 
the theory with interactions 



9{k) 



V{k, k', UJn) = D{k, k', Un) + 2 T(^n) ,3 



9{k') 



V^ik,k',ko)^Dj,{k,k',ko)- 



AB- {ko - itf 



a(fco) 



g(fcO 



(51a) 
(51b) 



For positive k^ we can reconstruct the retarded propaga- 
tor (j51bp from the temperature one (j51ap by making the 
substitution w„ —iko + t. Therefore, it is self-evident 
that the linear polarizabilty a{ko) can be extracted from 
the temperature transition matrix T(A;o). The rule is the 
following: 



a(fco) = -AT{-iko + e), 



for 



fc > 0. 



(52) 



Once we have determined the temperature transition ma- 
trix in the second order of perturbation ([55]) , we are able 
to find the polarizability in this order. For positive lo one 
gets 



AmA 



(4m2 - cj2) coth (^) - 4m [A(tj) + iT{Lo)] ' 



Functions A(w) and r(a;) are defined as real and imagi- 
nary part of the function h{Lj) = h{—iLU + e) 



A(w) = p 



dfc 



\k) 



fc2 



(53a) 
(53b) 



Finally, using the crossing symmetry of the propagators 
HH we get 



a(2)(a;) = 



AmA 



(4to2 - Lu^) coth - 4m [A{uj) + i sign(tj) T{lu)] ' 

(54) 

The last step is to use the correspondence formula 
to find a real physical atom polarizability. To do so one 
has to find a ratio of the partition functions Z/Z. In the 
lowest order of perturbation it is given by the formula 
(1112)) 



Z tanh(/3m) 
Z ^ 



tanh 



(^) 



tanh(/3m.) 
tanh ( 



(55) 



This all means that in the lowest order of perturbation 
theory the Fourier transform of the linear polarizability 
of a two-level atom has the form 



AmA-T 



Ara^ — up- — Am [At(cj) -|- i sign(a;) rT(ijj)] ' 

(56) 

where we introduced the effective temperature-dependent 
amplitude of the polarizability, as well as an effective shift 
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and width of the resonance 



Aj = tanh(/3m) 



^ _ tanh(/3m) | ^(2) 
tanh 



At(cj) = A{uj) tanh (^^^ , 
T,{lo) - r(c.) tanh f^) . 



(57a) 
(57b) 

(57c) 



It is very simple now to verify that in the low tempera- 
ture limit (/? ^ 00) we reconstruct the polarizability at 
zero temperature given by us before in pj]. Moreover, we 
find that our result predicts a small correction to the be- 
havior of the amplitude of the linear polarizability, well 
known from the theory of the magnetic resonance [35[. 
This correction is a second order outcome of the interac- 
tion and in the limit when g(k) — > we reconstruct the 
previous result. 

All physical properties of the linear polarizability of an 
atom at finite temperature directly follow from the equa- 
tion . As it is seen the magnitudes of the resonance 
shift and the width and the amplitude of the polariz- 
ability decrease when the temperature grows. It means 
that the capabilities of keeping the qubit under control 
at high temperature are extremely poor since the reac- 
tion of the system to the external electromagnetic field 
is weak (small amplitude), and the resonance is narrow 
(small resonance width). 



VIII. THERMAL LINE NARROWING 




FIG. 1: Temperature dependence of the relative shift and 
relative width of the resonance for two qubits with different 
energy gap. Energy gap Ai?'^' for the first qubit (grey line) 
is smaller than energy gap Ai?'^' for the second qubit (black 
line). Therefore, the slowing-down temperature is smaller in 
the first case than in the second case {Tc}^ < Tg^^). 



Qubit realization 


AE = 2m 


Ts 


Electron in quantum dot 


~ 50meV 


~ 1500 K 


^Si/a splitting in ^^^Hg+ 


167AieV 


~4.8K 


^Si/a splitting in "3Cd+ 


63/ieV 


~ 1.8K 


splitting in *''Rb+ 


28AieV 


~ 0.8K 



FIG. 2: Comparison of the slowing-down temperature Ts « 
2.5 X (temperature when life-time of exiting state becomes 
one order of magnitude larger than in zero temperature) for 
difi'erent physical realizations of qubits j3^, .38 . .39. . Ad . [IJ . 



Notice, that temperature dependence of the relative 
shift At/ A and the relative width Ft/F of the resonance 
in the lowest order of perturbation do not depend on the 
function g(k) but they depend only on the energy gap be- 
tween the qubit states AE = 2m (Fig. 1). It means that 
it is an universal property of physical qubits - it does 
not depend on the details of experimental qubit realiza- 
tions encoded in the coupling function g{k). Moreover, 
the interaction between qubits and the quantized elec- 
tromagnetic field is very similar to the interaction with 
the quantized vibrations of the crystal lattice (phonons). 
Therefore, one can imagine that the dependence on tem- 
perature in such a case would be very similar. This all 
means that relative temperature dependence of the shift 
and width of the resonance is universal for all qubits with 
freezed spatial degrees of freedom. 

From the eq. (|57cp one can find the sensitivity of the 
life-time of the exited state (the inverse of F) to the tem- 
perature for different experimental realizations of qubits. 
We define the slowing-down temperature Ts - the tem- 
perature when the life-time of the exited state becomes 
one order of magnitude larger than at zero temperature. 
The comparision of the slowing-down temperature for dif- 
ferent realizations is given in Fig. 2. 



In conclusion, let us concentrate on the problem of the 
width of the resonance. As was shown above, the reso- 
nance width unexpectedly decreases when the tempera- 
ture grows. Such a behavior is contrary to the natural ex- 
pectations. In typical situations we expect broadening of 
the resonance rather than narrowing. The best known ex- 
amples where the resonance broads with temperature are 
obviously heated gases where, according to the Doppler 
effect and the chaotic motion of the molecules, all spec- 
tral lines become wider when the temperature grows. 

In the literature there exist some theoretical predic- 
tions on the narrowing of the resonance known as the 
motional narrowing phenomenon. First theoretical anal- 
ysis of the motional narrowing was given in the funda- 
mental paper on nuclear resonance by Bloembergen et al. 

. The narrowing of the resonance is explained due to 
the influence of the thermal motion on the strength of 
the spin-spin interactions. It turns out that the thermal 
motion effectively leads to decreasing of the damping ef- 
fects induced by such interactions. The general problem 
of the motional narrowing was deeply studied and gener- 
alized in the language of the stochastic processes under 
the assumption that the system is subjected to a totally 
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chaotic fluctuating environment [43|- There exist lot of 

Ssical systems where the motional narrowing appears 

The origin of the narrowing of the resonance pre- 
dicted in presented paper is completely different. Notice 
that in the studied situation there is only one freezed 
qubit. Therefore there is no place for any influence of 
the Doppler effect or of any interaction with other qubits. 
This suggest that it should be called motionless narrow- 
ing rather than motional one, since it is clear that freez- 
ing of spatial degrees of freedom is a main source of the 
narrowing of the resonance. 



Acknowledgments 

I thank professor Iwo Bialynicki-Birula for his incisive 
questions, useful comments, and fruitful discussion. This 
work was supported by a grant from the Polish Ministry 
of Science and Higher Education. 



The argument leading to the formula (fT6|) is also valid for 
the systems with any A. Since the interaction Hamilto- 
nian Hi has the needed property (in the qubit subspace 
it is represented by a traceless matrix proportional to the 
Pauli ffx matrix), one has the important relation 



(A5) 



This immediately leads us to the nontrivial equality of 
the derivatives 



dZx 



dX dX ■ 
Therefore one gets a crucial relation 

Z — Zn = Z — Zn. 



(A6) 



(A7) 

It means that the difference between the partition func- 
tions Z and Z is the same as a difference between the 
partition functions for noninteracting systems Zq and Zq, 
which can be calculated directly. 



APPENDIX A: RELATIONS BETWEEN 
PARTITION FUNCTIONS 

In this appendix we derive a direct relations between 
the partition functions Z and Z. These relations are 
exact. They hold in any order of perturbation theory. 



1. Difference of partition functions 

Let us consider a family of the systems parametrized 
by A described by the following statistical Hamiltonian 



= ^0 + XTLi 



(Ala) 



and assume that each of those systems is in the quantum 
state described by a density matrix 



2a - tr 



(Alb) 



It is clear that for A = 1 we recover the situation de- 
scribed by the formulas ^ and for A = by the for- 
mulas ([6]). Let us also introduce the second-quantized 
counterpart of this family 



^A = A^o + XHi, 



Za = tr [e-^'^^ 



(A2) 
(A3) 



One can easily get a derivative of the partition function 
with respect to A for these two families 



dZx 
dX 

dZx 
dX 



-(3tT [e-^'^'-ni] 



Hi 



X 

P. 



Zx (XHi) 



(A4a) 



= ~^Zx {{XHi))x. (A4b) 



2. Ratio of partition functions 

From the general rules of the thermal quantum field 
theory we know that the partition function Z can be 
expressed as a sum of all vacum-like Feynman diagrams 



Zo (1 + <5) , 



where 



(A8a) 



(A8b) 



With the equation (jA7p one can find that 

Z ^ Za+ Z - Zo ^ Zo + ZoS (A9) 

and therefore the ratio of the partition functions is given 
by 

Z _Zo 1 + S 
Z ~ 



Zo 1 



Zo 



(AlO) 



This result is correct in any order of perturbation the- 
ory and it is easy to verify that the ratio Zo /Zo of the 
partition functions for noniteracting systems reads 



Zo 
Zo 



tanh(/3r7i) 



-f3m 



p/3m 



tanh 



(^) 



(All) 



To find an appropriate formula in a given order of pertur- 
bation one should expand the relation (|A10|) in powers of 
the coupling constant (the formfactor g{k)) and truncate 
the series at the appropriate position. For example, in 
the lowest (second) order one finds 



Z 
Z 



Zo 
Zo 



Zo 
Zo 



5(2) 



(A12) 



where S^^^ is represented only by one Feynman diagram 
(IB6a[) and it is given by (|B9p evaluated in the Appendix 



B. 
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APPENDIX B: SECOND ORDER 
CALCULATIONS 

1. Photon self-energy function 

The temperature self-energy function of the photon in 
the lowest order of perturbation theory is represented by 
the following Feynman diagram 



.p(2)(. 




(Bla) 



n' 

(Bib) 

Depending on the quantum statistics (see eq. dM])), the 
frequencies w„ and lo^' are 



2n7r 



{2n' + 1)tt 



(Blc) 



One can easily transform the relation (jBlaP to a simpler 
form 



_j oo 



n'— — oo 



1 



1 



iuJn' + iijJn + Tn iuJn' — m 

1 1 



IUJ„ 



iuj„ — m iujni + m 



(B2) 



The summation over n' can be done with the exploitation 
of the following mathematical identity 



E— — 

n ' 



-b 



^tanh(f ) -tanh(^) 
2 b-a 



(B3) 

which is valid for ujn' = (2n' + l)7r//3. 

The temperature self-energy function of the photon in 
the lowest order of perturbation theory has the form 



tanh + tanh f (rn + icj„) 



iuJn + 2m 



tanh {^] + tanh (^|(m - iw„) 
iun — 2to 



(B4) 



Conclusively, according to the periodicity of the hyper- 
bolic function and the relation (|Blc|) . the temperature 
self-energy function has the form 



4m 



■ tanh 



f3m 



(B5) 



2. Correction to the partition function 



The first correction to the partition function of the in- 
teracting system is represented by the following Feynman 
diagram 



^2 EE 







n n' 



dk dk' D{k,k',UJ.a) 



{v{k)S{uJr,, + iu^)V{k')SM} . (B6a) 



Summations run over the boson-like ujn and the fermion- 
like uj„' 



2m: 



LOrt 



i2n' + 1)tt 



(B6b) 



However, the summation over n' can be easily expressed 
by the photon self-energy function in the second order of 
perturbation theory (jBlaP evaluated before 



1 r°° p(-^ 



fc2 



4m / /3m 
= ^tanh — 



E 



dkg\k) 
1 



1 



+ Am? Lol, + /c^ 



(B7) 



The summation over n one can do similarly as in formula 
(|B2|) . For the bosonic frequencies (a;„ = 2n7r//3) the 
appropriate summation rule reads 



E 



lljJn — a lU), 



- b 



coth 



(f)-coth(f) 

b — a 



a) 

After performing the appropriate summations one gets 
the final expression for the first correction to the partition 
function 



tanh 



/3m 



oo 

dk 



k{k^ — Am?) 
pk 



2m coth ( — ) — k coth (/3m) 



(B9) 



At zero temperature limit (/3 oo), this correction is 
finite and negative. It is a simple exercise to find that in 
this limit it has a form 



p^oo Jq k[2m + k) 



(BIO) 



It is worthwhile to notice that even near the resonance, 
when k « 2m, the integral (jB9p is well defined, because 
the dangerous pole in the denominator is canceled by an 
appropriate term in the numerator of the expression. 
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